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HOLDER REGULARITY OF THE DENSITIES FOR THE 
NAVIER-STOKES EQUATIONS WITH NOISE 


MARCO ROMITO 

Abstract. We prove that the densities of the finite dimensional projections of 
weak solufions of fhe Navier-Sfokes equafions driven by Gaussian noise are 
bounded and Holder continuous, fhus improving fhe resulfs of Debussche and 
Romifo [DR14]. 

The proof is based on analyfical esfimafes on a condifioned Fokker-Planck 
equafion solved by fhe densify, fhaf has a "non-local" ferm fhaf fakes info ac- 
counf fhe influence of fhe resf of fhe infinife dimensional d 5 mamics over fhe 
finite subspace under observation. 


1. Introduction 

In this paper we improve the results of [DR14] for the law of the solutions 
of the Navier-Stokes equations with Gaussian noise in dimension three. We 
consider the problem 


( 1 . 1 ) 


u — vAu + (u ■ V)u + Vp = fi, 
V ■ u == 0, 


on the torus with periodic boundary conditions and driven by a Gaussian noise 
i]. In the equations above u is the velocity, p the pressure and v the viscosity of 
an incompressible fluid. It is known that the above problem admits global weak 
solutions, as well as unique local strong solutions, as in the deterministic case. 
Nevertheless the presence of noise allows to prove additional properties, such 
as continuous dependence on initial data [DPD03, DO06, FR06, FR07, FR08], as 
well as convergence to equilibrium [Oda07, Rom08]. See also the recent surveys 
[FR08, Debl3] for a general introduction to the problem. 


Our interest in the existence of densities stems from a series of mathematical 
motivations. The first and foremost is the investigation of the regularity prop¬ 
erties of solutions of the Navier-Stokes equations. We study here the regularity 
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properties of densities associated to the probabilistic distribution of the solu¬ 
tion, as existence and regularity of densities can be seen as a different type of 
regularity 

Our results concerns the existence of densities for suitable finite dimensional 
projections of the solutions of (1.1), and one reason for this is that in infinite 
dimension there is no standard reference measure (as is the Lebesgue measure 
in finite dimension), any choice should be necessary tailored to the problem at 
hand, and in our case we do not know enough of the problem. 

An interesting difficulty in proving regularity of the densities emerges as a 
by-product of the more general and fundamental problem of proving unique¬ 
ness and regularity of solutions of the Navier-Stokes equations. Indeed, a clas¬ 
sical tool is the Malliavin calculus, and it is easy to be convinced that it is not 
available here. Indeed, formally, the equation satisfied by the Malliavin deriv¬ 
ative is the linearisation of (1.1) and thus, estimates on the linearized equation 
are as good for the density as for uniqueness. We remark that in the case of the 
two dimensional Navier-Stokes equation, existence and smoothness of densi¬ 
ties for the finite dimensional projections of the solutions are proved in [MP06] 
with Malliavin calculus. 

This settles the need of methods to prove existence and regularity of the den¬ 
sity that do not rely on this calculus, as done in [DR14]. For other works in this 
direction, see for instance [DM11, BC14, KHT12, F1KF1Y13, HKHY14]. 

Existence of densities and their regularity in Besov spaces has been proved 
in [DR14] (see also [Roml3, Roml4c, Roml4a]), by extending and generalising 
a one dimensional idea from [FPIO]. Time regularity of the density has been 
proved in [Roml4b]. The method introduced in [DR14] is simple but effective 
and has been already used in other problems (see for instance [DF13, FoulS, 
SSSlSa, SSSlSb]) to prove the existence of densities. 

The results of [DR14] ensure that the density of the projection of solution at 
some fixed time on some finite dimensional sub-space is in the Besov space 
B|j^. Roughly speaking, this says that densities have (almost) integrable deriv¬ 
ative (see Section 2.1 for a short introduction to Besov spaces). 

In this paper we show a proof of Holder regularity of densities of finite di¬ 
mensional projections that is completely analytic and, unlike [DR14], does not 
rely on probabilistic ideas. In fact we follow a classical approach to existence 
and regularity of densities, namely the Fokker-Planck equation. The Fokker- 
Planck equation describes the evolution of the density of the ltd process solu¬ 
tion of a stochastic equation. Here we only look at a partial information on 
the solution (namely, a finite dimensional projection), thus we derive in Sec¬ 
tion 3 a Fokker-Planck equation with a "non-local" term that takes into account 
the effect of the dynamics outside the finite-dimensional space under observa¬ 
tion. The non-local term is indeed a conditional expectation and its regularity 
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is known only in terms of the unknown density itself. This makes our Fokker- 
Planck equation slightly non-standard. We re-derive in this framework the re¬ 
sults of [DR14] (see Proposition 4.1), we then prove the core result of the paper, 
namely boundedness of the densities, in Proposition 4.4, and finally the Holder 
regularity. 

Our proof of boundedness requires that, at least at the level of the Galerkin 
approximations we work with, we already know that the densities are bounded, 
possibly with bounds depending on the approximation (and so useless for the 
limiting problem). We derive these bounds on the approximations in the ap¬ 
pendix by standard results for hypo-elliptic diffusions and to do so we need to 
assume that the noise is "sufficiently non-degenerate" (see Section 2.3). With 
periodic boundary conditions the problem has been already thoroughly ana¬ 
lyzed in [Rom04] and this is the reason we mainly focus on the problem on the 
torus. There is in principle no limitation for the problem with Dirichlet bound¬ 
ary conditions, once smoothness of the densities at the level of approximations 
is settled. 

We believe that this is just a technical requirement (indeed it is not needed 
in [DR14]) that depends on the approach we have followed. Moreover, there 
is an inherent limitation in the Fokker-Planck approach, in that it is less flexi¬ 
ble fhan the method developed in [DR14] and thus caimot be used, in general, 
to evaluate the density of quantities that do not have an associated evolution 
equation, as for instance in [SSSlSa, SSSlSb], as well as for nonlinear functions 
of the solution of a diffusion process. The development of a probabilistic proof 
of fhe resulfs of fhis paper via a generalization of [DR14] is currenfly the subject 
of an on-going work. 


2. Main result 


2.1. Notations. We shall use the following nofations. If K is an Hilberf space, 
and F c K a subspace, we denofe by Ttp : K —)■ K the orthogonal projection 
of K onto F, and by span[xi,..., x^] the subspace of K sparmed by its elements 
Xi,..., Xn- Given a linear operator Q : K —)■ K', we denote by Q* its adjoint. 

We recall the definition of Besov spaces. Given f : —» R, define 


j=0 

and, for s>0, l^p^oo, l^q<co. 


(A}^f)(x) f(x + h) -f(x), 

n 

(A;;f)(x) := AUA;;-‘f)(x) = fL + iM. 

i —n 
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and for q = co. 


[fli 


sup 

IH^l 




LP 


where n is any integer strictly larger than s (the above semi-norms are inde¬ 
pendent of the choice of n, as long as n > s). Given s>0, l^p^oo and 
1 ^ q ^ oo, define 


B® 

p-q 


:= {f: 


LP 


+ < oo}. 


This is a Banach space when endowed with the norm ll^llB^q := l|f||LP + [fla^q- 

When in particular p = q = oo and s E (0,1), the Besov space B^^(R'^) 
coincides with the Holder space C^(R‘^), and in that case we will denote by 
II ■ II c| and ['Ici the corresponding norm and semi-norm. Notice that a more 
general definition of Besov spaces, that includes also the case s ^ 0, is based on 
the Littlewood-Paley decomposition. We refer to [Tri83, Tri92] for more details 
on the definitions we have given and the cormection with the general definition. 

For a Hilbert space K and a finite dimensional sub-space F of K, we shall 
denote by H(F), C^(F) and B® ^(F) the Lebesgue, Holder and Besov spaces, 
respectively, on F, when F is identified with a Euclidean space. 


2.2. The Navier-Stokes equations. Let H be the standard space of periodic 
square summable divergence free vector fields, defined as the closure of pe¬ 
riodic divergence free smooth vector fields with zero spatial mean, with iimer 
product (•, •)h and norm || ■ ||h- Define likewise V as the closure of the same 
space of functions with respect to the H^ norm. Let T7 l be the Leray projector, 
and let A = —TTl A be the Stokes operator. 

In view of Section 2.3 we specify in details an orthonormal basis of H. Let 
\ {0} and consider for every k G Z^ an orthonormal basis x^, of the 
subspace k-*- of R^ orthogonal to the vector k. Choose moreover xl,xl so that 
= xj^, i = 1,2. An orthonormal basis of H is given by £ = (e]^)t=i, 2 ,kez 3 , 
where ej^ = xj^ Clearly £ is a basis of eigenvectors of the Stokes operator 
A, and different choices of (xj^]i=i_ 2 ,kez 3 yield different bases. 

Define the bi-linear operator B : Vx V —)■ V' as B(u, v) = TTl (u ■ Vv], u, v G V, 
and recall that (ui, B(u 2 ,U 3 )) = — (U 3 , B(u 2 ,ui)). We will use the shorthand 
B(u) for B(u, u). We refer to Temam [Tem95] for a detailed account of all the 
above definitions. 

Let S : H —H be a Hilbert-Schmidt operator, then with the above notations, 
we can consider the following Navier-Stokes equations 

(2.1) du + (vAu + B (u)) dt = S dW, 

with initial condition u(0) = x G H, where W is a cylindrical Wiener process 
(see [DPZ92] for further details) in H. It is well-known [Fla08] that for every 
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X G H there exist^jrartngale solution of this equation, that^a filtered proba¬ 
bility space (n, P, a cylindrical Wiener process W and a process u 

wJA trajectories in C([0, oo); D(A)') fl i-i“(,([0, oo), H) fl oo); V) adapted to 

such that the above equation is satisfied with W replacing W. 


2.2.1. Galerkin approximations. Given an integer N ^ 1, denote by Hn the sub¬ 
space Him = span[ej^ : 1 = 1,2, |k| ^ N] and denote by ttn = tth^, the projection 
onto Him • It is standard (see for instance [Fla08]) to verify that the problem 

(2.2) du^ + (vAu^ + (u^)) dt = ttnS dW, 

where = TtiMBjTTiM-)/ admits a unique strong solution u’^ for every initial 

condition x’^ G H^. The proof is a based on the simple fact that in finite dimen¬ 
sion all norms are equivalent, thus given any finite dimensional sub-space F of 
H, there is Ci >0 such that 

IIttfAxIIh < Ci||x||h, ||7rFB(xi,X2 )||h ^ Ci||xi||h||x2||h- 

If X G H, x*^ = ttimX and Pj^ is the distribution of the solution of the problem 
above with initial condition x^ , then any limit point of (PN) iM^i is a solution of 
the martingale problem associated to (2.1) with initial condition x. In the rest 
of the paper we will consider only solutions of (2.1) of fhis type, as specified by 
the following definition. 

Definition 2.1. A solution of (2.1) with initial condition x G H is any process 
u with u(0j = X and with trajectories in C([0, ooj; D(A)') fl ([0, oo], H) fl 
co); V), such that there are a sequence of integers Ni< d oo and a se¬ 
quence Xic G Him,^ such that x-^ u(0), and u is a limit point, in distributions, of 

where u*^ is solution of (2.2). 


2.3. Assumptions. Consider again the Hilbert-Schmidt operator S : H ^ H. 
We will assume 


(2.3) 


S is diagonal in a basis £, 


where bases £ have been defined in Section 2.2. Under this assumption the 
operators A, S commute, £ is also a orthonormal basis of eigenvecfors of the 
covariance SS*, and the noise SW is homogeneous in space. This assumption is 
taken for the sake of simplicity (See Remark 2.2 below), in view of applying the 
results of [Rom04]. 

We shall also need the following global non-degeneracy condition, 

> let A = {k : (Sej^, ej^) > 0, i = 1,2}, then X is an algebraic system 
of generators of the group 


This assumption has been introduced in [Rom04] to ensure that Galerkin ap¬ 
proximations (2.2), with N large enough, are hypo-elliptic diffusions. We will 
use this fact to deduce, see Theorem A.l, that the solution of (2.2) has a smooth 
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density with respect to the Lebesgue measure on Hn . Notice that the assump¬ 
tion essentially requires that the noise is stirring all directions in H, albeit indi¬ 
rectly We believe this assumption is technical and depends on the way we have 
proved our results, namely to ensure that the computation in the next section 
are rigorous. As already mentioned in the introduction, a probabilistic proof of 
our main results would definitely get rid of this assumption. 

We will consider densities for the projections of solutions of (2.1) over finite 
dimensional sub-spaces of H. Given a finite dimensional subspace F of H, we 
consider the following conditions. 


(2.5) 


F is the span of a finite subset of £, 


where the basis £ is the same of assumption (2.3), and that the noise is non¬ 
degenerate on F, namely. 


( 2 . 6 ) 


7rpSS*7rp is a non-singular matrix. 


This condition has been introduced in [DR14, Roml4b], and amounts to say that 
the covariance has full range in F or, in different words, that the noise is directly 
stirring all directions of F. 


Remark 2.2. In general, there is nothing special with the bases £ provided by the 
eigenvectors of the Stokes operator and our results would work when applied 
to Galerkin approximations generated by any (smooth enough) orthonormal 
basis of H. On the other hand we are using the results of [Rom04] and the 
setting with the bases £ is the most suitable. 

One could rather work with a general basis of eigenvectors and assume that 
the spectral Galerkin approximations are hypo-elliptic diffusions, thus ensur¬ 
ing the conclusions of Theorems A.l and A.4. We have preferred to proceed 
with the explicit version of the assumptions. Similar considerations apply for 
the problem on a bounded domain with Dirichlet boundary conditions. 


Remark 2.3. The two non-degeneracy assumptions given above are, in a way, 
independent. An easy example where assumption (2.4) holds while assumption 
(2.6) does not is provided in [RXll], where all but the "low modes" are forced 
by the noise and (2.4) holds. But if F contains low modes components, (2.6) is 
not true. 

On the other hand, fix ko G and set 


S = <(•> + ^n(-, 4ko)enko> 

rL=l 

SO that the set X introduced in assumption (2.4) is X = {nko : n G Z, n 7 ^ 0}. It 
is easy to check that, formally, the Navier-Stokes dynamics (2.1) is closed in the 
subspace span[ej^ : i = I, 2, k G X], and in particular is not hypo-elliptic. If F is 
sparmed by elements of X, then (2.6) holds. 
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We do believe (and this is the subject of a work in progress) that also in this 
case the projections on F have FFolder densities. Notice that Proposition 4.1 
below only ensures that these densities are in B j 

2.4. Main result. The main result of the paper is as follows. 

Theorem 2.4. Consider equation (2.1) and assume S satisfies assumptions (2.3) and 
(2.4). Let ube a solution of (2.1) as in Definition 2.1 and let F be a finite dimensional 
subspace ofH satisfying conditions (2.5) and (2.6). For every t > 0 denote by fF(t] the 
density o/Ttpujt) with respect to the Lebesgue measure on F. Then for every T > 0 and 
every oc G (0,1), 

SUpt^t‘'+«’||fF(t)||c«<00, 

(0,T) 

where d = dim F. 

In particular the density fp is bounded and FFolder continuous for every ex- 
ponenf a. < 1. The proof of fhis fheorem will be given af the end of Section 4. 


3. Formulation of the conditioned Fokker-Planck equation 


Under our assumptions (see Theorem A.l) we know that, if u’^ is a solution 
of (2.2), fhen for every t > 0 fhe law of u*^ (t) has a smooth density Fn (tj with 
respect to the Lebesgue measure in the Schwartz space. Then it is a standard 
fact that fN satisfies a Fokker-Planck equation, fhaf in our nofafions reads 

(3.1) 9tfN = ^^nFn + Vn ■ ((vAnx + BN(x])fN), 

where 

^/Ng = Tr(SNSND^g), 
and Vn' is fhe divergence on Hn. 

Fix a subspace F of Hn such fhaf conditions (2.5) and (2.6) hold and consider 
the projection Tipu’^ of u*^ on F. The marginal densify is given by 


fF,N(t,x') 


Fn (t, x' + x") dx", 

F^N 


where F-*-^ is the space orthogonal to F in Hn . In the sequel we will understand 
X G Hn as X = (x',x") with x' G F and x" G F-'-'^. We wish to derive now an 
equation satisfied by fF,N- By infegrafmg the equation (3.1) over F-’-^, it is not 
difficult to see that 

(3-2) Stfp.N = + Vp ■ (Gp^Nfp.N), 

where sfy — Ttp^e/NTtp (the resulting operator is independent of N due to the 
assumption (2.3) that the covariance is diagonal in the Galerkin basis), and 


(3.3) Gf^nIFjXO = VTtpAx' + EjTTpBN (u*^ (tjjlTTpu'^ (t) =x']. 
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Indeed, set for brevity Gn (x) = vAnx + Bn (x), then 


Vn - (GnFn) — 

Vf - (ttfGnFn) + 

Jf^n 

F-^N 


F^N 


VpJ-N • (tTp^nGnIn) 


where the second integral in the displayed formula above is zero by integration 
by parts, and the first integral can be reinterpreted as 


Vf ■ (ttfCnIn) dx" — 


F-Ln 


F-Ln 


Vf ■ (7tFGNfF-LN|F(f,^l’^')fF,N(t,x')) dx' 


= Vf - (E[7rFGN(u(t))|7rFu(t) = x']fF,N(t,x')). 


Here is the conditional density of u^(t) given 7rFu’^(t], and Vf- 

is the divergence on F. An additional simplification, due to the fact that An is 
also diagonal in the Galerkin basis, yields (3.3). Similar but simpler computa¬ 
tions show that the contribution of .c^nIn/ when averaged over F-'-'^, is the term 

Lemma 3.1. Let N he large enough (that F c Hn) and u*^ he a solution of (2.2). Then 
Gf,n (t) G LP (F; fF,N dx')/or every p > 1 and t > 0. Moreover, for every T > 0 there 
is Cz — C 2 (p, v,T,F, ||u’^(0)||h) > 0 such that 


(3.4) 


Sp(T)^ := sup 

[ 0 ,T] 


|GF,N(t,x')l^fF,N(t,x') dx' < C2. 


Proof Fix p ^ 1 and t > 0. Let fpi-NiFli^A'I^') be the conditional density of 
u’^ (t) given Ttpu’^ (t), then by the Holder inequality. 


|GF,N(tA 0 r ^ 


IttfGn (x', x"]|'"fF±N|F(t,x"|x') dx", 


F-Ln 


hence 


|GF,N(t,x')pfF,N(t,x') dx' < 


|7rFGN(x]|’"fN(t,x] dx 




HJ 


= E[||7rFGN(u^(t))|| 

^ C3(vPE[||u^] + E[||7rFBN(u^]) 

^C4(vPE[||u^(t]||^]+E[||u^(t)fJ’]) 

^ C2, 

where we have used the fact that on F all norms are equivalent and the (uniform 
in N) estimafe (A.l). □ 

Remark 3.2 (exponential bound for Gf,n)- A slightly (although useless so far) 
better estimate can be obtained using (A.2). By the Jensen inequality 






F-^N 
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hence, as in the proof of the lemma above, and by (A.2), 

gA|GF,N(t,x')l dx' ^ ^ E[e^^5(i+I|u^(t)||^)]^ 

Jf 

that is finite for A small enough. 

Remark 3.3. The equation (3.2) can be recast in a way that shows more explic¬ 
itly how the contribution of modes in F-’-^ enfer in the evolution of f f,n • More 
precisely, for x e Hn, 

7rFBN(x) = 7rFBN(x') -F 7rFBN(x',x") -F ttfBn (x", x') -F7rFBN(x"), 


and so 

Gf,N (f, = vAfx' + Bf(x') + E[7rFBN (tTfu'^ (t), Ttpi^(f)) I TTfU*^ (t) = x'] 

-F E[7rFBN (TtpXN (f) , TTfU*^ (t) ) I TTfu'^ (t) = x'] 

-F E[7rFBN (7tp±N u'^ (t)) I ttfu’^ (t) = x'], 

thaf is fF,N solves an equation analogous fo (3.1) wifh an additional term that 
takes into account the influence of the evolution of modes from . 


3.1. Integral formulation of (3.2). Due fo our assumptions (2.5) and (2.6) on F, 
the operator .s/f is elliptic and with constant coefficients. Let p*" be its kernel, 
then the integral formulation of equation (3.2) is given by. 


(3.5) fF,N(t,x') = 
where Xq = Ttfu’^ (0). 


rt 

0 


(VfpL, 


-k 


(Gf,N (s, •)fF,N (s, •]) (^0 ds, 


Proposition 3.4. Under assumptions (2.5) and (2.6) on ¥, for every p with 1 ^ p ^ oo 
and every n ^ 1 there is Cs > 0 such that for every h G F with |h| ^ 1, and t > 0, 

(3.6) ||p[||lp ^ C6t“^ IIVPpIIlp ^ 

(3.7) ||AK)3,f||„ $C5t-A(lAE)'', IIA^Vp'IIl, 5:cst-lA+l)(lA.^)’', 

where q is the Holder conjugate exponent ofp and d = dim F. 


Proof The operator £/y is a second order elliptic constant coefficients operator, 
so its kernel Pp is, up to an invertible linear change of variables, the standard 
heat kernel. In particular, by parabolic scaling, p[(x) = t“‘^/^pp(x/\/t). With 
these observation at hand the proof of the proposition follows from similar 
compulations for the heat kernel. The computations are elementary and thus 
omitted. □ 
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4. Holder regularity 


In this section we prove our main Theorem 2.4. Prior to this, we give a 
proof of the result of [DR14] using the Fokker-Planck formulation (see Propo¬ 
sition 4.1). Notice that, consistently with [DR14], we do not need assumption 
(2.4) to do so. Then, as an intermediate step, we prove boundedness of the den¬ 
sities. This is, not surprisingly in fact, the crucial step and with boundedness at 
hand the Holder regularity follows easily. In order to obtain L°° bounds we need 
to know though that densities are smooth (albeit with bounds not necessarily 
uniform in N). We will take care of this in the appendix. 

Notice finally that, unlike in [DR14], we have not been able to derive better 
bounds for stationary solutions. This may have a twofold reason. On the one 
hand the higher summability needed to derive better bounds for stationary so¬ 
lutions are not as good as those of Lemma 3.1. On the other hand the method of 
[DR14] exploits non-trivially correlations in time that are harder to use in this 
framework. 

4.1. Basic Besov regularity. In this section we wish to give a proof of the Besov 
regularity of the density fp alternative to [DR14], using the formulation with the 
conditioned Fokker-Planck (3.5). This is to emphasize that in [DR14] Lemma 3.1 
is (implicitly) used only with p = 1 and thus that there is space for improve¬ 
ment. This will be then the subject of the rest of the section. 

Proposition 4.1. Let f be a subspace of H satisfying (2.5) and (2.6). Then there is 
C/ > 0 such that for every N ^ 1, Yq G Hn, and t > 0 



where fp^N is the density of the solution of (2.2) with initial condition Xq. 

Proof Clearly HIp^nIIli = 1/ so we need to estimate only the Besov semi-norm 
Leth G F, with |h|p ^ 1, then by Proposition 3.4 



We have by (3.5) that 


^ (Gp,N(s)fp,N(s)) ds, 
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where Si is defined in (3.4). Therefore, 



1A 


\/t — s 



«A=(i+Si(t))ihi, 

Vl At 


and, by definition, [fF,N]B;„ <4c6(lAt) 


□ 


In the above proof we have only used that Gf n G L^(F). As long as we 
only know this, the previous result is essentially optimal. Indeed, the term 
Vf • (GF,NfF,N] is, roughly, in hence, by convolution with the heat ker¬ 

nel, we can expect that fF,N is at most in (that is very close to what we have 
proved). 

Since by Lemma 3.1 the quantity Si is uniformly bounded in N, in the limit 
N —)■ oo we can derive a result for solutions of the infinite dimensional problem, 
as in [DR14]. 

Corollary 4.2. Given a weak martingale solution u of the Navier-Stokes equations as 
in Definition 2.1 and a finite dimensional sub-space f of H such that conditions (2.5) 
and (2.6) hold, for every t > Othe random variable 7tFu(t] admits a density in 
with respect to the Lebesgue measure on F. 

Proof Let u be a solution of (2.1) according to Definition 2.1. Then u is a limit 
point, in distribution, of the sequence (u’^)n^i of solutions of (2.2). By Propo¬ 
sition 4.1 and the embedding of B( ^ in LT’ (for a p > 1 depending on the di¬ 
mension of F), the densities fF,N (f) of ttfu*^ (t) are uniformly integrable. We can 
thus find sub-sequences, that we will keep denoting by (u’^ )n^i and (fF,N )n^i/ 
such that u*^ converges in distribution to u and Ff.n (t) converges weakly in L^ 
to the density fF(t) of 7tFu(t). The Bj ^ bound on If follows now easily because 
of the weak convergence. □ 

4.2. Boundedness of the densities. As an intermediate step in the proof of our 
main result we prove that under our assumptions the densities are bounded. 
This is the crucial step and having boundedness at hand, Flolder regularity then 
is not difficult to prove. 

Given a finite dimensional sub-space F of H satisfying conditions (2.5) and 
(2.6), an integer N ^ 1 large enough that F c Hn, and a solution u*^ of (2.2), 
define for every T > 0 and cx > 0, 


3"f,n(T]:= sup t‘^||fF,N(t)||oo 


(4.1) 


te[o,T] 


Lemma 4.3. Under assumptions (2.3) and (2.4), ifN^l is large enough (that F c 
Hn), there is ocq ^ f such that ifu^ is a solution of (2.2) and fF,N (■) is the density of 
ttfU*^ (■)/ then (T) < oofor every T > 0. 
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Proof. Fix k > N. By Theorem A.4 we know that there is cci^ > 0 such that 
M := supt x(l A t)‘^'‘(l + |x|’^)fisi(x] < oo. Hence 


(1 At)‘^'"fF,N(t,x') = (1 At)‘^'‘ fN(t,x) dx" A 

J F-'-N 


< 


pJ-N 1 + |X 


//Ik 


(1 A t)‘^'^(l + |x|^)fN (x) dx" ^ CgM. 


Therefore Tp< oo. □ 

The next result provides a uniform (in N) estimate of Tp ^(T). 

Proposition 4.4. Under assumptions (2.3) and (2.4), given T > 0 and M > 0, ifV is a 
subspace of H satisfying (2.5) and (2.6), there is Cw = Cwiy, T, F, M.) > 0 such that if 
N is large enough (that F c Hn), and ifu^ is a solution of (2.2) with Hu*^ (0)||h ^ M., 
then 

3^f%'(T) ^ CIO. 

Proof Fix T > 0 and a solution u*^ of (2.2) with initial condition Xq, and set 
Xq = tt^Xq. By (3.2), for every x' G F and t e (0, T], 


rt 


(4.2) 


|fF,N(t,x')l ^ llPtlll 


||Vpp_3 -k (Gf^n (s))||l° 


ds. 


Fix e G (0, and let d = dim F. On the one hand, by Proposition 3.4, 


pt—te 


rt—te 


||Vpp_3 ★ (Gf,N (s)fF,N (s] 


ds A 


IIVpL 


(4.3) 


sIlL” 
pt—te 


^ Sim 

^ Cii(te 


'GF,NfF,N)(s)||Li ds 
C6 


0 


/4- 

(t-s) 2 


ds 


,-Ad-i) 


Sim. 


On the other hand, if cx > 0, p G (1, and q is such that ^ ^ = 1, by the 

Holder inequality and again Proposition 3.4, 

(4.4) 


ft 

||VPp_s ★ (Gf,N (s]fF,N (s) 


ds ^ 


t—te 


imPt-sIkp |KGF,NfF,N)(s)||Lq ds 

C6 


t—te 


rt 


^3^F%m^Sq(T) 


t-te SP (t — S) 2 


ds 




1 d 1 _ a d 


Ci2e2 2qt2 p 2qTF,N(T]pgq(T). 
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Take |3 > -I V + p ~ then by using together (4.3) and (4.4) into (4.2) we 
obtain 

Tp%(T] + 

+ Ci2e^-^TP+^-t-^Tp«Nm^9q(T]. 

We can now prove that Tp (T) < oo. Indeed, let (Xq be the exponent given by 
Lemma 4.3, so that Tp “pj (T) < oo. If ao = f there is nothing to prove. Otherwise 
we can take cci = fV(^ + -^ — t) and, by (4.5), Tp (T) < oo as well. By the 
choice of p (that gives q > d), it turns out that ai < ocq. It is also clear that by 
iterating the above procedure with exponents 0,2, (X3, ... , in a finite number of 
steps we will obtain the exponent ^ and thus that Tp (T) < 00 . 

Finally, we prove the uniform bound on Tp ^(T) < 00 . Consider again (4.5) 
with a = |3 = -I, and use the Young inequality, 

3^f,n(T) ^ C6 + Cne-^'‘'-'^\/Tgi(T) + Ci2e^-^x/Tgq(T)Tp%"(T)E 

^ C6 + Cne-^‘‘i-'^\/Tgi(T) + l(ci2e^-^\/Tgq(T])^ + -Tp^iJlT). 

q V M / p r,iM 

Since p > 1 we deduce that ^ C(T, SilT), Sq (T)] and the conclusion of 

the theorem finally follows. □ 

The singularity in time of the L°° norm in the previous result clearly originates 
only from the singularity in the initial condition. It is reasonable then that, when 
we look for an initial distribution for u*^ with a smoother law, the same result 
should hold with a smaller power for the time. As an example, we give a direct 
proof of the L°° bound of the density of the (unique, see [Rom04]) invariant 
measure. Denote by kp n its density, then as in Section 3, the density satisfies 
the equation 

-.e/pkF^isi + Vp ■ (Gp^Nkp^N) = 0. 

Notice that when the averaging in (3.3) is done with respect to the invariant 
measure, Gp,n does not depend on time. Lemma 3.1 holds though and Gp,n G 
LP (F; kp^N dx') for every p ^ 1. 

Proposition 4.5. Under assumptions (2.3) and (2.4), given a subspace F 0 / H satis¬ 
fying (2.5) and (2.6), and N ^ 1 large enough, let kp n be the density of the unique 
invariant measure of (2.2). Then there is C 13 > 0 depending only on v, F and some 
polynomial moment of the invariant measure in H such that ||kp^iM ||loo ^ Cu. 

Proof Let d = dim? and denote by gp the Green function of jYj=. As in Propo¬ 
sition 3.4, gp can be obtained by an invertible linear transformation from the 
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Poisson kernel. Hence Vgp ^ c^lxl^ We have 
kF,N (ic') = —(Vgp) T*r (GF,NkF,N)(K-0 

— ~(^9FlBe(x')) * (GF,NkF,N)(K-0 — (VgFlB|(x')) * (GF,NkF,N ) (K-0 

= [!] + §) 

where e > 0 will be chosen later. Fix q > d and let p be its Holder conjugate 
exponent. Then 

IQI ^ ||VgF||LP(Be( 0 )||GF,NkF,N||Lq ^ CisCP S q || kF,N || l°° % 

and likewise 

Q,-%(r 1— — 

101 ^ ||VgF||LP(B|(0)||GF,NkF,N||LP ^ 9pll^F,N ■ 

where Sp/ Sq are the quantities in (3.4) but computed on the stationary solution 
(so there is no need to evaluate the supremum in time). 

Collecting the two estimates above and choosing e = ||kF,N yields 

|1<-F,N (^01 ^ CiseP ^^Sq ||kF,N “I-Sp ||kF,N 

^ Ci7(Sq + Sp) ||kF,N ||oo^ , 

and hence the statement of the proposition, since ||kF,N ||l°° is non-zero and fi¬ 
nite by Lemma 4.3. □ 

Clearly both results above immediately extend to solutions of the infinite di¬ 
mensional problem that are limit point of Galerkin approximations as in Corol¬ 
lary 4.2. 


4.3. Holder regularity. Using the boundedness of the densities proved so far 
we can give an estimate on Gf,n different from (3.4). Indeed, if p > 1 and t > 0, 


|GF,N(t)fF,N(t)||LP 


|GF,N(t,x')fF,N(t,x')P dx' 


(4.6) 


^ iifF,N(t)ir ^ 


I 

rd/2, 


|GF,N(t,x')|PfF,N(t,x') dx' 


^ Sp(t)Tp;N(t)H 


where q is the conjugate Holder exponent of p. In particular it is elementary to 
verify that Gp.Nfr.N ^ l-’^([0, T]; (F)) for every T > 0 and every p, r ^ 1, with 

I -F p > d. In other words Gp.Nfr.N has classical summability. 


Proposition 4.6. Under assumptions (2.3) and (2.4), given a subspace ¥ of H satisfy¬ 
ing (2.5) and (2.6), cx e (0,1), T >0 and M > 0, there is Cjg = Cjg(v, a, T, F, M.) > 0 
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such that if'N is large enough (that F c F-In), and ifu^ is a solution of (2.2) with 
llu’^ (0)||h ^ M, then 

sup II c« ^ Ci8, 

te[o,T] 

where d = dimF and (■) is the density ofu^ 


Proof. Fix a G (0,1), T >0 and a solution of (2.2) with initial condition Xq, 
and set Xq = TtpXo. By (3.2), for every x' G F, t G (0, T] and h G F with |h| ^ 1, 


II AHfF,N (t) ||l°o ^ ||AHPtl|l 


ft 


IKAnVpt-s^ ★ (GF,N(s)fF,N(s))||L“ ds. 


Fix e G (0, |] and set q = On the one hand, by Proposition 3.4, 


•t—te 


II (AhVp[_s) A (Gf,N (s)fF,N (s] 


ds ^ 


•t—te 


iiahVp; 


^9i(T)|hr 


sl|L“ 
't—te 

0 


[GF,NfF,N)(s)||i_i ds 
C6 


fj. 3 

(t-S) 2 


ds 


^Ci9(te)-^gi(T]|hr 


On the other hand, by the FI older inequality, the estimate (4.6) and again Propo¬ 
sition 3.4, 


II (AhVPt-s) * (GF,NfF,N ||l“ ds ^ 


t—te 


||AHVp^_5||Lp||GF,NfF,N||Lq ds ^ 

t—te 


pt 


C6 


d . d I 1 

t-te S2p (t — Sj 


1 A 


ds A 


\/t — S' 

^ C2it-^'i-i’jp%'(T)?gq(T)|h|'-A 

where p is the conjugate Fiolder exponent of q. Since by definition 1 — = a, 

the two estimates above together yield 

||AHfF,NW||L» ^ + Ci9(te)-^t‘'+«-i5gi(T)|hr + 

and therefore 

sup(o(‘i+-nfF,N]Bs») ^ C22(1 + \/Tgi(T] + x/rg a 
[0,T] ’ 


The conclusion follows since ^ = Cg^. □ 

Proof of Theorem 2.4. Let u be a solution of (2.1) according to Definition 2.1. Then 
u is a limit point, in distribution, of the sequence (u’^ )n^i of solutions of (2.2), 
and the density Ff of ttfU is a limit point for the weak convergence in of 
(fF,N)N^i by Corollary 4.2. 
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Proposition 4.6 above states that the sequence (fF,N)N^i is bounded in C^(F], 
hence by the Ascoli-Arzela theorem each sub-sequence admits a further sub¬ 
sequence converging uniformly on compact sets of F (and weakly in since it 
is also uniformly integrable) to a C^(F) function. This proves the theorem. □ 


Appendix A. Technical estimates 


In this appendix we derive some (non-uniform in N) quantitative bounds on 
the density Fn of u*^ that are necessary as an intermediate step in the proof of 
our main result. To this end we recall that the solutions of the Galerkin approx¬ 
imations (2.2) satisfy 




f-T 

(A.l) E 

sup ||u^(t)||H^ + V 

IIu^HvIIu^IIh'’'' dt 


40,T] J 

0 


^ C23(p,^,T, ||u^(0)||h,S], 


for every N ^ 1, u’^ (0) G FI, T >0 and p ^ 1, and the number on the right hand 
side depends on N only through Hu*^ (0)||h- In fact there is a stronger estimate: 
there is A > 0 such that for every IM ^ 1, u’^ (0) G FI and T > 0, 


(A.2) E 

exp (A sup 1 u’^l ^ 4- vA 

1 u’^ 1 / ds) 


L ^ [0,T] 

0 


^ C24(y,T, ||u^(0)||h,S), 


See for instance [FlaOS] for details. 

The main point in our estimates of Section 4 is that, in order to prove uniform 
(in N) bounds in L°° for the marginal density fp.N/ we already need to know 
that fp N is regular. This in principle does not follow immediately from the 
regularity of (see Example A.3 below) unless the joint density Fn is in the 
Schwartz space, as shown in the result below. 


Theorem A.l. Under assumptions (2.3) and (2.4), /or every xq G FIn the density Fn 
of the solution u’^ of (2.2), with initial condition xq, is in the Schwartz space. 

Proof The conclusion follows immediately from [Nua06, Proposition 2.1.5] once 
we show that u*^ (t) is a non-degenerate random variable (see [Nua06, Defini¬ 
tion 2.1.1]), namely that u’^ (t) G D°°, where is the space of random variable 
with Malliavin derivatives of arbitrary order, and that for 

every p ^ 1, where MiM^t is the Malliavin matrix of u’^ (t). 

The fact that u*^ (t) G D°° for every t > 0 can be proved as in [Nua06, Theorem 
2.2.2]. The quoted theorem requires that the coefficients should have bounded 
derivatives of any order larger or equal than one. The coefficients of our equa¬ 
tion (2.2) are second order polynomials, and it is not difficult to verify, by go¬ 
ing through the proof of the quoted theorem, that the boundedness condition 
can be replaced by (A.l) to prove existence and uniqueness, and by (A.2) to 
prove existence of the Malliavin derivatives, since the equations for the Malli- 
avin derivatives are linear. 
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Finally, the fact that the inverse of the Malliavin matrix has all polynomial 
moments finite follows as in [Nua06, Theorem 2.3.3] (where again we use (A.2) 
to replace the boundedness condition on the derivatives), once the Hormander 
condition (see [Nua06, Section 2.3.2]) holds. Under assumptions (2.3) and (2.4) 
the Hormander condition follows from [Rom04, Lemma 4.2]. □ 


From the previous theorem we immediately deduce the following conse¬ 
quence, whose proof is based on elementary computations. Notice that, in view 
of Example A.3, the fact that the density fN is in the Schwartz space is crucial 
for the proof of boundedness of fF,N • 

Corollary A.2. Under the same assumptions of the previous theorem, the marginal 
density n is the Schwartz space. 

Example A.3. It is easy to construct a density 4> £ C°°(R‘^) fi such that at 

least one of its marginals is not bounded. Indeed, take d = 2 and let cp G C°°(R^) 
such that cp ^ 0, Supp cp C (-i, ||(p||li = Land ||(p(0, OUli + ||(p(-,0)||li > 0. 

Set for every k e Z^, 

A, ( ^ i ki ^ 0, ka ^ 0, 

^ |0, otherwise, 

and (j) = X!k ~ Then 4> is a smooth, bounded density (up to re¬ 
normalization), but both marginals are unbounded. 

We give a more quantitative version of the previous theorem by giving the 
explicit dependence of the supremum norm of the density with respect to time. 
This is necessary for the proof of Proposition 4.4. Similar bounds can be also 
given for the derivatives of the densities, although we do not need these esti¬ 
mates in the article. 


Theorem A.4. Under assumptions (2.3) and (2.4), let Xq G Hn and let fN be the 
density of the solution of (2.2) with initial condition xq. Then for every k ^ 1 there 
is (Xi< > 0 such that for every T > 0, 

sup ((1 A IxI’^jfN (x)) < oo. 

te(o,T),xeHN 


Proof By the previous theorem we know that f n is in the Schwartz space, so we 
only have to understand the singularity at t = 0. By [Nua06, Proposition 2.1.5] 
we have the following formulas for the density. 


Dn 

fN(t,x] = e|^]^ .(t)>Xj}HN,t 
j=l 


[n! ^NNj(t)>xd^NN,i(t)<Xi}TlN,t 
iA 
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for X G Hn and t > 0, where Dm = diniHiM, Hn t G D°° is a suitable random 
variable arising from integration by parts, and ._Dn the compo¬ 

nents of (t) with respect to a basis of Hn . By the Cauchy-Schwartz inequal¬ 
ity, the representation formulas above and (A.l), 

(1 -|- |x|’^)fN(x) ^ C28lE[(l + (t]||H)|Ht|] ^ 

< C29E[(1 + ^ c30E[|HtP]A 

for every k ^ 1. By [Nua06, formula (2.32)], given p ^ 1 there are |3,y > 1 and 
integers n, m ^ 1 such that 

where Du*^ (t) and MN,t are respectively the Malliavin derivative and the Malli- 
avin matrix of and || ■ ||k,p, tor k ^ 1 and p ^ 1 is the norm of the 

Malliavin-Sobolev space see [Nua06, Section 1.2], By [Nua06, Theorem 

2.2.2] , together with our exponential bound (A.2), as in the proof of the previ¬ 
ous theorem, we see that 

sup ||2)u^(t)||DN,T < OO, 

[0,T] 

and that the singularity at t = 0 originates, not unexpectedly, from the inverse 
Malliavin matrix, or equivalently, from the smallest eigenvalue AN,t of Mnu- 
Indeed, | det ^ hence 

IldetM^yirPio) 

To conclude the proof, we briefly go through the proof of [Nua06, Theorem 

2.3.3] and point out only, for brevity, the arguments where the dependence in 
time arises. To show the tiniteness of the moments of A]i]\, we look for the tail 
probabilities P[AN,t ^ e] at 0. By [Nua06, Lemma 2.3.1] 

E[^N,t ^ e] ^ C 32 e^^^'^ sup P[v'''MN,tV ^ e] -|-P[|MN,tl ^ 

|v|=l 

The second term on the right hand side is not the source of singularity in time 
and can be easily bounded since MN,t has all the polynomial moments finite. 
The first term is bounded using the Norris lemma ([Nua06, Lemma 2.3.2]). In¬ 
deed for every q ^ 2 there is C 33 > 0 such that 

sup P[v^MN,tV ^ e] ^ £ 336 *^, 

|v|=l 

if e ^ eo, for a suitable Eq > 0. The value of Cq is identified through the afore¬ 
mentioned Norris lemma, and here one has to look for the singularity in time. 
It is elementary to see that, just by going through the proof of the lemma, that 
eo ~ (1 A t)*^, for some a that depends on q. By this it follows that 

E[A-P'"^] ^C34(lAt)-‘^P°A 
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as needed. 


□ 
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